Critical and tricritical exponents of the Gross-Neveu 
model in the large- N { limit 

Hiroaki Sugiyama* 
Department of Physics, Tokyo Metropolitan University, 
(N ■ Minami-Osawa, Hachioji, Tokyo 192-0397, Japan 

H- S 



Abstract 

The critical and the tricritical exponents of the Gross-Neveu model are calculated 
■^J- \ in the large- iVf limit. Our results indicate that these exponents are given by the 

C*") ' mean- field values. 
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1 Introduction 



The thermodynamical properties of QCD is relatively unexplored subject which is nev- 
ertheless important in its own right []]]. Furthermore, it is currently receiving renewed 
interests because of the heavy ion collision experiments started at Relativistic Heavy Ion 
^ Collider (RHIC) which is expected to have potential of probing into transitions to the 

quark-gluon plasma phase. 

It is, however, highly nontrivial problem how to extract the thermodynamic quanti- 
ties of QCD from these experiments. The critical exponents, for example, are defined in 
equilibrium statistical mechanics and it is not obvious how they can be revealed in such 
energetic collision experiments. Therefore, it may be of help if we have a model which has 
the similar thermodynamic properties as QCD and at the same time simple enough to be 
tractable analytically. 

It was shown in their original work that the Gross-Neveu model has the same physical 
properties as QCD, asymptotic freedom, dimensional transmutation, and the spontaneous 
chiral symmetry breaking 0. The thermodynamics of the Gross-Neveu model in the large- 
Ni limit is discussed by Wolff who worked out the structure of the phase diagram of chiral 
symmetry breaking on chemical potential-temperature plane, and in particular found a 
tricritical point ||. Furthermore, Chodos et al. have argued recently that by adding a 
pairing interactions the extended model admits a superconducting phase and possesses a 
very similar phase diagram with that of two-flavor massless QCD with chiral symmetry 
breaking and Cooper pairing phases M. Therefore, the model is a good candidate for the 
theoretical laboratory of QCD for studies its thermodynamic properties. 
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In this paper, we study the thermodynamics of the original Gross-Neveu model and 
calculate the critical and tricritical exponents. It can be regarded as a preparatory study 
toward the discussions of how to extract thermodynamic quantities from robable dynamical 
processes in nonequilibrium or moderately equilibrium environments. A work has already 
been initiated aiming at pursuing such a goal 0. 

In this paper, we restrict ourselves to the large-iVf limit of the Gross-Neveu model. Since 
it is a (1 + l)-dimensional field theory, nontrivial phase structure at nonzero temperature 
is only possible in the large- iVf limit due to the Mermin- Wagner theorem M. The nature 
of the phenomenon and its interpretation as an "almost spontaneous symmetry breaking" 
has been discussed by Witten [[FJ . We think it sufficiently illuminative of the fact that the 
Gross-Neveu model in the large- Nt limit provides a consistent theory which is usable as a 
theoretical laboratory for QCD. 

In Sec. 2 we review the Gross-Neveu model at finite chemical potential and temperature. 
In Sees. 3 and 4 we calculate the critical and the tricritical exponents of the Gross-Neveu 
model. The final section is devoted to conclusions. 

2 The Gross-Neveu model at finite density and tem- 
perature 

The Gross-Neveu model is a (1 + l)-dimensional model whose Lagrangian density at finite 
fermion density may be described by 

C GN [ip,ip;iJ l ,mo) = $i(i$ - mo)ipi + — — {ipiipif + ^aoipi (2.1) 

where /i denotes the chemical potential and mo the bare mass. The index i denotes flavor 
of fermions and runs from 1 to N{. Note that when we set mo = the model has an 
invariance under the discrete chiral transformation if) — > —^75, ip — > 75^. The theory is 
renormarizable, by virtue of that the model is (1 + l)-dimensional one. 

To make the model easy to handle we introduce an auxiliary field a by adding a term 




£ aux = -77T \G ~ m + -Tffpifpi (2-2) 



to the original Lagrangian density £gn- Considering the fact that JVaexp (i /e$r£ aU x) is 
an irrelevant constant we can use the Lagrangian density 

£ = C GN + £ aux = if>i(i$ + /ry - a)4>i - —(a - m ) 2 (2.3) 

without changing the dynamics of the original Lagrangian density. 

For uniform vacuum we define the effective potential V e s(a) as T c s = —iNf (Jd 2 x) V e s, 
where r e ff is the effective action. Now we take iVf infinity. In the limit we can construct 
V e Q from tree and one fermion loop diagrams as shown in Fig. |l|. We obtain 

V eS (a } /1, mo) = —(a - m ) 2 + i J — — In (-(fc + /i) 2 + k\ + a 2 ) (2.4) 
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up to a constant term. From the equation of motion for a we see 

cr c i = (a) = m - — (^i(x)if>i(x)) ■ 



(2.5) 



The vacuum expectation value a c \ is the physical mass of fermions and also order parameter 
for the discrete chiral symmetry breaking. 

So far we dealed with zero temperature case, but from now let us put the system into 
thermal bath of finite temperature T. That is achieved by replacing fc - m tegral with sum 
over Matsubara frequencies as follows: 



> > i>T, 
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(2n + l)mT 



(2.6) 



where the Boltzmann constant is set to be unity. After the replacements the effective 
potential becomes 

V eS (a, fjt, T, too) = 7rr(cr-m ) 2 



2A 



+ i J^ E ^ln[-{(2n + l)7uT + / i} 2 + ^ + a 2 



■ (2-7) 



Now we must do renormalization, because the second term of (|2.7| ) clearly diverges. Naively 
speaking, 1/A and rriQ are the parameters which we can use for renormalization. But it 
is better to regard 1/A and m /A as independent parameters. We define a renormalized 
coupling constant Ar as follows: 
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where a' denotes an arbitrary renormalization point and A is a momentum cutoff. In fact 
the divergence in (|2.7|) can be absorbed by only the coupling constant renormalization. 
Therefore, m /X is finite and needs not to be renormalized. Now it is convenient to define 
m F by 



TOp = cr' exp ( 1 — — ) 
V Ar/ 



(2.9) 



which represents the dynamically generated mass of fermions for ji — T — m® = 0. We 
obtain the gap equation 

dV eS (cr, fjt, T, m ) 
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(2.10) 



where 7^ denotes Euler's constant (~ 0.5772), see || for detailed calculation. The phase 
diagram on \i-T plane for mo = was obtained in || (Fig. |2|). In the region surrounded by 
the curve ABCO the discrete chiral symmetry is spontaneously broken. On the curve AB 
there occures a second-order phase transition, while BD is the first-order phase transition 
line. Thus, the Gross-Neveu model in the large- Ni limit has a tricritical point B. 



3 The critical exponents of the Gross-Neveu model 

We study the thermodynamical properties of the Gross-Neveu model in the large- Nf limit. 
At first we obtain the critical exponents analytically. In order to do this we set /i = for 
the time being. Then the gap equation ( |2.10| ) becomes 
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(3.1) 



In the last equality we have used high temperature expansion which can be used for the 
region 1/T < n/\cr c i\. We can use the expansion near the second-order phase transition 
point, because a c \ is almost zero and T is nonzero. The critical temperature T c (the 
temperature of the point A in Fig. ||D has been given [|10] by solving 
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(3.2) 



and the result is 



T r 



e 7E rvj 0.566933 mp. 
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When the temperature is T c and mo goes to zero, the gap equation becomes 

dV cS 







da 



CT = CT cl 
(J = 
T=T C 



m 7C(3) 3 
A 8tt 3 T 2 c1 



-I^T" 4 



(3.3) 



(3.4) 



and the behavior of the order parameter a c \ is given by 

0"cl 



1/3 1/5 

-^ m = m • 



(3.5) 



Therefore, we obtain a critical exponent 5 = 3. On the other hand, when mo = and the 
temperature closes to T c from bellow, the gap equation becomes 
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cff 
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and the behavior of a c \ is 
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(3.6) 



(3.7) 



Thus, we obtain a critical exponent (3 = 1/2. There exist other critical exponents a and 7 
defined by 
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which govern the behavior of the specific heat Cy and the susceptibility % near the critical 
point. Because we have already obtained two exponents 5 = 3 and (3 = 1/2, we can obtain 
them by using the scaling relations 



a = 2 -(3(8 + 1) = 0, j = /3(S-l) 



(3.9) 



4 The tricritical exponents of the Gross-Neveu model 

From now on we reinstall the chemical potential \i and study the tricritical exponents of 
the Gross-Neveu model in the large- Nf limit. See, for example, [|llj for an extensive review 
of the tricritical phenomena. We solve the gap equation ( |2.10| ) numerically and see the 
behavior of the solution a c \ near the tricritical point. First of all we need to know the 
tricritical point /i tc and T tc . They are obtained by numerically solving the equations 
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d*V eS 



da* 



cy OO 

~ £Re 



J^2 



tr=mo=0 ■* n=0 

and the results are given [|| by 



{(2nH-l)7T + ^/T} c 



^tc ~ 0.608221 m F , T tc ~ 0.318328 m P . 



(4.2) 



(4.3) 



We can achieve high temperature expantion of fl2.10|) [p|j which can be used for the region 
1/T < 7r/(|cr cl | + |/i|). Although the condition is satisfied near the tricritical point, the 
expantion does not enable us to calculate the tricritical exponents analytically because the 



power series of 1/T is not of a for fi ^ 0. Therefore, we calculate the tricritical exponents 
numerically. 

At the tricritical point the behavior of a c \ will be charactarized by a tricritical exponent 
S tc as a c \ ~ m Q tc for sufficiently small m . Our numerical computation indicates 



1 m da, 



cl 



Stc cr c \ dm, 



A» = Mtc 
T =T tc 
mo 2i 



0.20 (4.4) 



as shown in Fig. §. We roughly estimate the error of 1/S tc to be less than 0.01. 

Next we set mo = and investigate the behavior of a c \ near the tricritical point. Though 
there are many paths to approach to the tricritical point and be two exponents which relate 
with the paths, at first we choose simply the path of T = T tc . As approach to the tricritical 
point along the path, the order parameter will behave as a c \ ~ (/i c i — //) * ! with a tricritical 
exponent $. Our numerical computation indicates 



0, 



•i ^ 



Htc - \x da, 



cl 



tfci d(fj, tc - n) 



M - Mtc 
T =T tc 
run = 



0.25 (4.5) 



and is shown in Fig. |[ We roughly estimate the error of (3i to be less than 0.01. 

There should be another exponent corresponding other paths, but to obtain it is rather 
difficult. For example the path of \i = /i tc gives just same exponent as fy and in fact any 
linear paths do also. Here remember that the curve AB in Fig. ^| causes the second-order 
phase transition and that the curve AB and BC are parts of the solution of an equation 
( f4.1|) . The path of T = T tc which crosses the curve ABC had given an tricritical exponent 
f3i ~ 0.25, therefore the path along BC (in broken phase) may give another exponent which 
must exist. For the reason we define the parameter s which parametrizes the curve BC 
in units of rrip so that s = at the tricritical point B. Then, the order parameter a c \ 
behaves as a c \ ~ s^ with a tricritical exponent (3 t near the tricritical point. Our numerical 
computation indicates 



„ s da c 
Pt 



r-^/ 



ct c i ds 



0.5 (4.6) 



as shown in Fig. [|. That is indeed the other exponent than f3i that we have looked for. But 
the accuracy of j3i is less than that for /3 t because of rather intricate way of approaching to 
the tricritical point. Though there are more tricritical exponents a\ and at for the specific 
heat, 7; and jt f° r the susceptibility, we can obtain them from the scaling relations 



aj = 2-A(<*tc + 1) - 1.5, 


li = Pi(Stc - 


-1)-1, 


a t = 2-(3 t (5 tc + l)~-l, 


It = Pt(Stc - 


-1) ~2 



(4.7) 

It is remarkable that a t has a negative value. It indicates that the specific heat is continuous 
at the tricritical point as a function of s, which is different from what is expected in usual 
second-order phase transition. 



Finally, let us confirm the validity of our numerical calculations for the tricritical expo- 
nents. By essentially the same computer program for the tricritical exponents, we numeri- 
cally calculate the critical exponents 5 and (3 which have been already obtained analytically. 
The results are shown in Fig. |] and Fig. [7[ and these results agree with analytical ones very 
well. Therefore, we believe that our numerical calculations for the tricritical exponents are 
also accurate enough. 

5 Conclusion 

In this paper, we have calculated the critical exponents of the Gross-Neveu model analy- 
tially and the tricritical exponents numerically in the large- N{ limit. The results we have 
obtained are as follows: 

the critical exponents : 5 = 3, j3 = — 
the tricritical exponents : — = 0.20, (3i = 0.25, j3 t = 0.5. 

The accuracy of our numerical calculations of the tricritical exponents are checked against 
the values of the critical exponents which we can obtain analytically. The errors in the 
tricritical exponents are thereby estimated as less than a few %. 

Our results indicate that all the critical and tricritical exponents of the Gross-Neveu 



model in the large- N{ limit are given by the mean-field values [0. The reason may be 
that the fluctuations around a mean field are suppressed in the leading order of 1/iVf 
expansion. Unfortunately, to calculate the next order of 1/iVf does not help because the 
critical temperature becomes exactly zero at finite Nf in the Gross-Neveu model; it is 
believed that the result of applies to any field theories in (1 + l)-dimensions with short 
range interactions. 

Then, what is the significance of the computation of critical exponents we have carried 
out in this paper? It is to prepare for the study of the question of how equilibrium 
thermodaynamical quantities can be signaled and extracted from the experiments, as we 
have discussed in Introduction. However, we have to keep in mind, despite their similarities, 
the following differences between the Gross-Neveu model and QCD. Namely, we have no 
reason to expect the same numerical values of the critical and the tricritical exponents, 
and furthermore their physical nature may not be the same. It is because the Gross-Neveu 
model in the large- Nt limit have the mean-field critical exponents, while QCD with two 



massless flavor is believed to belong to the universality class of 0(4) spin model [IB ]. We 
hope that these differences do not seriously disturbe the validity of the Gross-Neveu model 
as a theoretical laboratory for studying thermodynamic properties of QCD. 
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Figure 1: The diagrams which contribute to the effective potential in large- N{ limit. Only 
the terms up to fourth order in a are explicitely exhibited. 



Figure 2: The phase diagram of the Gross- Neveu model in the large- Nf limit and the forms 
of the effective potential each for tjiq = |1. The point B is the tricritical point. 
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Figure 3: The tricritical exponent l/5 tc is plotted as a function of mo/m-p. Here, /i and T 
are set to be // tc and T tc respectively. In the vanishing mo limit 1/S tc approaches to 0.2. 
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Figure 4: The tricritical exponent (3i is plotted as a function of (/i tc — /i)/mp. Here, m = 
and T is set to be T tc . As /i closes to /i tc , A approaches to 0.25. 
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Figure 5: The tricritical exponent j3t is plotted as a function of s which parametrizes the 
curve BC in Fig. ^| in units of rnp so that s = at the tricritical point 5. The unit of s is 
mp. Here, m is set to be zero. As s approaches to the tricritical point, (3 t approaches to 
0.5. 
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Figure 6: The critical exponent 1/5 is plotted numerically as a function of mo/m-F- Here, 
fi = and T is set to be T c . In the vanishing itlq limit 1/5 approaches to 0.33 which is 
contrasted with the analytical result 1/5 = 1/3. 
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Figure 7: The critical exponent (3 is plotted numerically as a function of {T c — T)/m-p. Here, 
m and \x are set to be zero. As T closes to T c , /3 approaches to 0.5 which is contrasted 
with the analytical result /3 = 1/2. 
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